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ERROR ANALYSIS OF A DIFFUSE INTERFACE METHOD FOR 
ELLIPTIC PROBLEMS WITH DIRICHLET BOUNDARY CONDITIONS 

MATTHIAS SCHLOTTBOM^ 


Abstract. We use a diffuse interface method for solving Poisson’s equation with a Dirichlet 
condition on an embedded curved interface. The resulting diffuse interface problem is iden¬ 
tified as a standard Dirichlet problem on approximating regular domains. We estimate the 
errors introduced by these domain perturbations, and prove convergence and convergence 
rates in the H^-norm, the L^-norm and the L°°-norm in terms of the width of the diffuse 
layer. For an efficient numerical solution we consider the finite element method for which 
another domain perturbation is introduced. These perturbed domains are polygonal and 
non-convex in general. We prove convergence and convergences rates in the H^-norm and 
the L^-norm in terms of the layer width and the mesh size. In particular, for the L^-norm 
estimates we present a problem adapted duality technique, which crucially makes use of the 
error estimates derived for the regularly perturbed domains. Our results are illustrated by 
numerical experiments, which also show that the derived estimates are sharp. 


Keywords: diffuse domain method, elliptic boundary value problems, embedded Dirichlet conditions, im¬ 
mersed interface, domain approximations, finite element method 
AMS Subject Classification: 35J20 65N30 65N85 

1. Introduction and main results 

This paper considers the approximate solution of the following model problem by a diffuse 
interface method: Find u G such that 

— Au = f in n \ r, ttir = fl'ir on F. (1) 

Here, the function / G models volume sources in a convex polygonal domain C M"", 

n = 2,3, and the function g G defines the values of u on an interface F G 

where F C is a closed manifold of co-dimension one, i.e. for n = 2 a curve, or a surface if 
n = 3. We assume that the interface F separates 12 into two domains Q = Di U F U D 2 , where 
F = dDi, see Figure [T| 

The analysis of ([^ is well-established. For instance, if g does not depend on u, Q can 
be separated into two independent Dirichlet problems on Di and D 2 respectively, and the 
theory for the Poisson equation with Dirichlet boundary conditions applies, cf. [26^ 129] , 
Alternatively, one may formulate Q as a Dirichlet problem on Q constrained by tt = 5 on F 
and u = 0 on 5D, which leads to a saddle-point formulation, see e.g. [niET]. 

The numerical approximation of Q has been investigated intensively. A first set of numeri¬ 
cal algorithms relies on a triangulation of D which is sufficiently aligned with the interface, i.e. 
F is approximated by a polygon, the segments of which are edges (or faces) of the elements; 
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Figure 1 . Sketch of the geometry. 




see e.g. laiiaEaiiii 1211111 ]. The construction of such a triangulation might be expensive 
or difficult in practice. Furthermore, if F = T{t) depends on time, problems similar to Q 
need to be solved in each time step, and the triangulation has to be updated accordingly. In 
addition, one has to interpolate the data on the varying meshes in this case. Therefore, a lot 
of research has been conducted to construct accurate methods employing meshes which are 
not aligned with F but possess a “simple” structure and are fixed throughout the simulation; 
see for instance the immersed boundary method [l 2 | , the immersed interface method ISlEZj, 
immersed finite elements |36[ l3^ HH] , the fictitious domain method laiiziiiHiiioi, the unfit¬ 
ted finite element method la tta EH [SI, the finite cell method m, unfitted discontinuous 
Galerkin methods [9|, or composite finite elements [.SOI I38j . and the references provided there. 

In this work we will focus on a diffuse interface method for solving Q, see for instance 
[IlliaEHlIlllEZlESlIlsl. In this method the sharp interface condition u = g onT is replaced 
by suitable conditions on tt — 5 on a diffuse layer centered around F. Similar techniques 
have also been applied for solving coupled bulk-surface differential equations [ 1 ] or surface 
differential equations dniEoi. 

The constraint tt = (7 on F is equivalent to the condition ||u — 5 ||x, 2 (p) = 0. In order to relax 
this condition on the sharp interface, we define the signed distance function 

|-|-dist(x, r), X G T>2, 

and we let 5 : M —)• M be such that S{t) = t for |t| < 1 and S{t) = sign(t) for |t| > 1 . 
The width of the diffuse layer is characterized by a positive parameter e. This induces a 
regularized indicator function of Di 


XDi{x) ^ Uj%x) = J(l + 


X G fl. 


2' ■ e 

Formally dF = |Vx£)|dx « |Va;^|dx. Using the e-tubular neighborhood = supp(|Vcu^|) 
of r, this leads to the following approximation for integrals along the interface 






The reader might find a more detailed derivation of this approximation and other choices 
of S in |17j . We further notice that, to make this approximation well-defined, we need the 
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function g to be defined on the diffuse layer. If g is defined on F only, one has to use a suitable 
extension; for instance a local extension is given by + (ir(a^)Vfir(3^)) = g{x) for a: G F, i.e. 
g is constant off the interface. Thus, replacing the sharp interface constraint by the diffuse 
interface constraint \u — g\‘^ dx = 0, which amounts to n = 5 on 5^, we are concerned with 
the following Dirichlet problem: Find G Hq(Q) such that 

—Au^ = f infl\5^, = g in 5^. (2) 

Note that the particular choice of io^ is not important as long as = supp(|Va;^|) is a 
e-tubular neighborhood of F, i.e. methods using double obstacle potentials to regularize the 
indicator function of XD whl essentially lead to the same method. 

The main purpose of this paper is to estimate the errors introduced by this diffuse interface 
method; (i) on the continuous level and (ii) in a finite dimensional setting when using the finite 
element method, see below. The first result in this direction is the following approximation 
result on the continuous level. 


Theorem 1.1. Let f G and g G and let u be the solution to Q, and let he 

the solution to Q. Then there exists a eonstant C > 0 independent of e such that 
1 1 

-||m — u^\\l^(q) H—^||Vt6 — Vu^\\l^(q) < (^(11/1^2(0) + \\9\\m{n)) ■ 
e ye 

The numerical approximation of Q is still not straight-forward as the (sufficiently exact) 
integration over is basically not easier than the integration along F. In order to obtain an 
efficient numerical scheme let 7/i be a shape regular triangulation of Ll, and let 

= U T (3) 

T£Th-S^nT^0 

denote the union of all elements having non-empty intersection with S^; see Figure below 
for one instance of and Sf^. Here, h = max{diam(r) : T G Th} denotes the mesh-size 
parameter. We are then concerned with the following approximate Dirichlet problem: Find 
^e,h g Hq{LI) such that 


= f mn\Sl 


= g 


= 9 in St. 


(4) 


Note that Q is still formulated on the continuous level. The difference to ([^ is that we have 
replaced the domain Ll\S^, which has the same smoothness as D \ F, by D \ 5^, which is a 
polygonal domain. Problems similar to Q for the approximation of Q have been considered 
previously e.g. in [12] for n = 2; see also nn for very general approximating domains. We 
will employ similar techniques as in m to derive Theorem |1.2| from Theorem o 

Theorem 1.2. Let f G L^(D) and g G and let u he the solution to ([^, and let 

be the weak solution to Q. Then for some C > 0 independent of e, h, and 6 there holds 

^ _l_ ^11“ “ ^^’^11^2(0) -I- --^==||Vn - < C'(||/||i2(Q) -I- ||5 '||f/ 2 (q)), 

where 6 = max{diam(T) : T n dS^ / 0, T G Th}- 


Concerning the numerical approximation of we denote by I 4 = PiiTh) C ^o(^) n 
standard finite element space of continuous piecewise affine functions. For the incorporation 
of the interface values we denote by L^g G Pi{Th) H the nodal interpolant of g, and let 

= {vh e Vh : VhiT = ih9)\T for T n cS*" 7 ^ 0}. (5) 
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The Galerkin approximation G ^hg defined by the variational problem 


Vvhdx= / fvhdx lox aXlvh^V^Q. 


(6) 


By construction Sf^ is aligned with Th- Therefore, integration of the bilinear form in Q can be 
performed by standard techniques without additional variational crimes. Note also, that the 


extension of g is only needed on . Using Cea’s lemma. Theorem 1.2 and problem adapted 
duality arguments, we will also prove the following theorem. 

Theorem 1.3. Let f G L‘^{Ll) and g G and let u he the solution to Q, and let 

G he defined via Q Then for some C > 0 independent of e, 6, k, and h there holds 

^||m - - Vu^||2^2(f2) < C'(||/||x,2(f2) + \\g\\H^(n)), 


where A = ^/e + 5 + ks + h, K = max{diam(T) : T n 
max{diam(T) : T G Th\- 


T G Th\ and h = 


In view of the above theorems, the diffuse interface method investigated here is robust 
in the choice of e and h. In particular, e might be chosen arbitrarily small such that the 
overall approximation error is governed by approximation errors from the finite element ap¬ 
proximation. In this case the described method is similar to the method described in [34j . 
However, we do not impose the restriction that our triangulation is <5-resolved or e-resolved 
in the sense of [Ml Definition 3.1], and, also, dSf^ C basically induces two polygonal inter¬ 
faces approximating T, whereas in [Mj, one polygonal interface is defined. The ^-estimate 
in Theorem 1.3 has a similar form as the estimate proven in |34l Theorem 4.1], but here the 
value of 6 is already determined by the mesh size. The value of 5 in |34j stems from the 
assumption that the triangulation is (5-resolved. If the hnite element grid is aligned with T, 
i.e. if the vertices defining dSf^ are points on dS^, stronger approximation error estimates can 
be derived [niiaiiiiiii]. 

Let us compare our results with those in the literature about diffuse domain methods for 
solving the Dirichlet problem Q in Di. In [T7j the following variational problem has been 
considered: 


Find G such that 


Vvf^-Vvoj^ dx+ — [ u^v\Vuj^ 
! « Jq 


dx = [ fvuj‘^dx + — f gv\Vui^ 
Jn OL Jq 


dx 


(7) 


for all V G , where is a weighted Sobolev space, which consists of functions in the 
weighted space having also weak derivatives in this space; details are provided 

in m- Note that functions in can develop singularities where = 0. The analysis 
in m requires a careful balance between a and e. Provided u G the optimal 

choice of a = gives a rate \\u% — = 0{£^^^) in [IT]. Let us point out, that, 

under additional assumptions, a rate 0(e) for the Lfi error could be obtained for n < 2. If 
u ^ W^’°°{Di) but only in which is essentially what we need here, the convergence 

rates of m get worse and depend on the dimension n due to embedding theorems. However, 
the computations in [33] show a convergence for the L^-error of order 0(e) for a similar diffuse 
domain method, where the diffuse interface condition is incorporated by a penalty method 
similar to Q; see Remark |5.4| below. In [25| an L°°-error estimate of order e^ for s < 1 


arbitrary has be proven for a diffuse domain method with a double well regularization for 

XDi 


and n = 1. As proven in Theorem 4.3 below, the diffuse interface method considered 
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here yields ||tt — = 0{s + 5). We consider only Dirichlet boundary conditions 

here, for Neumann or Robin boundary conditions let us refer to where under appropriate 
regularity conditions on the data even superlinear convergence has been shown for the diffuse 
domain method on the continuous level. 

The rest of this paper is structured as follows. In Section we introduce further notation 
and recall some solvability and regularity resu lts fo r 0 and 


1.2 


and L^-error estimates, which prove Theorem [lT} Theorem | 
completeness we also give an L°°-error estimate. Galerkin approximations are investigated 


In Section]^ we derive H 
Section]^ For 


IS proven m 


m 


Section where also Theorem |1.3| is proved. Our findings are supported by numerical 
examples in Section Section gives some conclusions. The paper closes with an appendix 
recalling some estimates for the errors introduced by the diffuse interface method. 


2. The Dirichlet problem revisited 


Let D C M”' be some bounded domain with Lipschitz boundary. By Lp‘{(X) we denote the 
Lebesgue space of square integrable functions. Accordingly, for an integer A: > 1, H^(p) = 
VF^’^(D) is the set of functions in L^(D) having also weak derivatives up to order k in L^(D). 
These spaces are endowed with the standard norms, i.e. 


Il2(Q) 


= / |u|^dx. 




k 

E 

| a |=0 




L2(Q)! 


where a G Nq is a multi-index. We will write Vtt for the gradient of u. Hq{Ci) is the closure 
in of infinitely often differentiable functions with compact support in Cl. We recall the 

Poincare inequality [21 6.26] 


Il2(Q) < 


diam(D) 


II^^IIl 2 (o) for all u G iLo(^)- 


Hence, ||Vu||j^ 2 (q) defines an equivalent norm on Hq{Ci). By we denote the restriction of 
a function u to D; if the context is clear, we will write u instead of u^q. 

In the whole manuscript we let eo > 0 be such that (i) the projection of the diffuse layer 
onto r is well-dehned, see (28), and (ii) (30) holds. Furthermore, we let 0 < e < Eq. We 
denote by C a generic constant which may change from line to line, and, in particular, C is 
independent of e, 5 and h and the data. Under the given assumptions, we have the following 
existence and regularity result, which follows from a combination of results from standard 
theory on elliptic equations 


Lemma 2.1. Let F G , f G L^(D), g G Then there exists a unique solution 

u G H^in) to © such that G H^{Di), z = 1,2, and 

+ II^|DiII_h2 (£,^) -I- \\u\£)2\\h^{D2) ^ C'(||/||2,2 (q) -I- \\g\\H'^(n))- 


For the solution of the diffuse interface problem we have 


Lemma 2.2. Let F G , f G L^(D), g G H‘^{Ll). Then there exists a unique solution 
G FAq (D) to ([^ such that \ 5^), i = 1, 2, and 

+ ll“fDi\5^llH'RDi\5‘^) + ll“fD2\5^lli^RR>2\5‘^) - C'(||/||l2(q) -|- 11511/12(0)). 
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For later reference, we note that the solutions u to 0 and to satisfy 

— u = 0 on r. 


( 8 ) 


3. Interface problems involving 

We start with error estimates for the solutions to the model problem Q and its diffuse 
interface reformulation Q. First, we derive estimates in the Fl^-seminorm, a nd t hen we 
derive a corresponding L -error estimate via the Aubin-Nitsche lemma. Theorem 0 is then 
a consequence of Theorem 3.1 and Theorem 3.2 below. Since Q is a linear equation, we may 
assume without loss of generality that 5 = 0 in this section. For the sake of simplicity, we 
let D = Di and D‘^ = D\S^ in the rest of this section, and perform the analysis for this 
domain. In particular F = dD. The remaining estimates on D 2 can be derived analogously. 
Furthermore, in slight abuse of notation, we denote by tt G Hq{D) r\H‘^{D) the restriction of 
the solution to Q to D. By G n we denote the corresponding restriction 

of the solution to Q. In particular, there holds 

— Au = / in D, u = 0 on dD, (9) 

—Au^ = f in n = 0 on dD’^. ( 10 ) 


Hence, we have for all v G Hq{D) that 


/ Vu-'Vvdx= / fvdx. 
Id Jd 


( 11 ) 


Extending = 0 on F^ = L) \ and integration by parts shows that for any v G Hq (D) 


Vu‘^-Vvdx= / Vu^-Vvdx= / fvdx+ dnu'^vda. 


( 12 ) 

Jd JD^ J D^ JdD^ 

Here, n denotes the exterior unit normal field to dD^, and dnU^ denotes the normal derivative 
of Subtracting © and ([T^, we obtain 


/ V{u —u‘^)-Vvdx = / fvdx— dnu'^vda for all u G Fig (I?). 
Id Jr, JdD'^ 


(13) 


Theorem 3.1. Let f G L‘^{D), and let u G Hq{D) n H^{D) be the solution to Q, and let 
G Hq{D^) n be the solution to (10). Then 


\\Vu-Vu^\\L2^D)<CVi\\f\\LHD)- 

Proof. Using (1^, we see that u — ^ Hq{D) is a valid test function for 


i.e. 


||Vu — Vu^||^ 2 f£)') = / f{u — u^)dx— / dnU^{u — u^) da. 

^ Jr^ JdD^ 


( 14 ) 


Using Theorem A.l and ([^, we obtain 

[ /(u-M^)dx < Ce||/||i2(r^)||Vu-Vu^||i2(^). 
JFe: 

Furthermore, an application of the Cauchy-Schwarz inequality yields 


IdD^ 


dnU^iu - U^) dcr < \\dnU^\\L^{dD<^)\W “ 


L'^(dD^)- 
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The term \\dnU^\\L'^[dD^) can be estimated in terms of ||/||l 2 (£)) using Lemma A.3 and Lemma 2.2 

and Lemma A.2, which gives 


For the other term, we employ 

11^^ - < C^/e\\Vu - V«^||i2(£,). 


Using these estimates in (14) completes the proof. 


Using the Aubin-Nitsche lemma |21) . we obtain 

11 “ - ^^IIl 2 (d) ^ - ^u^Wl^d) inf llVz - Vr;||i 2 (^), 


where ^ G Hq{D) n H^{D) denotes the solution to 

—Az = u — in D, 

Using V = e Hq{D‘^) n defined as the solution to 

-Az^ = u — u‘^ in D‘^, 


□ 


(15) 


which exists by Lemma 


2.1 


in (15), we obtain as a direct consequence of Theorem 3.1 the following statement 


Theorem 3.2. Let the assumptions of Theorem \3. 1\ hold true. Then 

Ik “ < C£||/IIl2(£)). 

4. Interface problem involving Sf^ 

As a hrst step towards a practical numerical scheme, we let 7^ be a shape regular trian¬ 
gulation of n, cf. 13 Dehnition 4.4.13]. Replacing the domain Ll\S^ by 11 \ 5^, where Sf^ 
is dehned in ([^, yields an important change in the geometry. Firstly, the distance of points 
in dSl to F is not constant anymore, whereas dS^ = {\dr\ = e}, and secondly, for fixed Th, 
dSl is only polygonal. Therefore, we cannot rely on Fl^-regularity of the dual problem to (Q. 
As in the previous section, we consider only the case g = D = Di and = Di\Sf^ in 
detail; the remaining cases follow with similar arguments. Similar to the previous section, we 
denote by G the solution to (10) with replaced by where 

6 = max{diam(r) : rn5^ / 0,r G Th}- We consider (|^ in weak form, i.e. let G 77g(L>|) 
be the solution to 


/ Vu^’^-Vudx= / fvdx for all u G Fig (H^). 

Id? Jd? 


(16) 


Since D? is a non-convex polygonal domain in general, the solution to (16) is not regular 


enough to repeat the arguments of the previous section. Our error estimate is based on the 
following observation 


Ds+S c Di C C D. 


(17) 


Theorem 4.1. Let f G L^{D), and let G Hq{D) n H^{D) he the solution to (10), and let 
g L[Q{Dff) be the solution to (16). Then there holds 


\\Vu^ -Vu^’^h^iD) <CV6 \\fh 2 ^n) 

with 6 = max{diam(r) : T 0 dS^ / 0, T G %_}■ 
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Proof. In view of we have that HQ{Dff) is closed in where functions in 

are extended by zero. Using Cea’s Lemma [15] and G we obtain 

||Vu^ - < inf ||Vu^ - < ||V«^ - < CVSWfh^^j,^, 


where we also used Theorem 3.1 with D replaced by and u replaced by u^. 
An application of the Aubin-Nitsche lemma [21] yields 
M - u‘ 


l 2 (D) < C\\Vu^ - inf - Vv\\l 2 (^d). 


where G Hq{D^) H is the unique solution to 


□ 


(18) 


-Az^’^ = — u' 


€,h 


in 


The first term on the right-hand side of (18) can be estimated using Theorem 4.1 The second 
term can be estimated using a similar reasoning as in the proof of Theorem 4.1 Summarizing, 
we have the following statement. 


Theorem 4.2. Let the assumptions of Theorem \4.1\ hold true. Then there holds 

\\u^ — < C'(5||/||2,2(£)) 

with 6 as in Theorem\4. 1\ 


Theorem 1.2 is now a direct consequence of Theorem 4.1 and Theorem 4.2 in combination 
with Theorem 3.1 and Theorem |3.2[ respectively. Complementing the result in |25|, we also 
give an L°°-error estimate; see 


Theorem 4.3. Let the assumptions of Theorem 4-1 hold true, but let additionally f G LP{D) 
for some p > n. Then there holds 

\\u^ - U^’^\\l‘^(DI) < C''5||/IIlp(D) 


with 6 as in Theorem 4-1 For the solution u G Hq{D) n of Q, there holds 

ll^e _ < C{e + (5)||/||lp(z)). 

Proof. We proceed as in |33|. We observe that w = u — u^'^ G is the weak solution to 

—Aw = 0 in Df, w = u on dDf. 


Using the weak maximum principle |261 Theorem 8.1], we see that 

II £ h\\ 11 11 

\\u-U’ \\l'=°{DI) < \m\L°°{dDl)- 

In view of pS] ) below, every x G dDf^ can be written uniquely as x = x -|- di'{x)n{x), where 
X = pr{x) G dD and n{x) = V(ir(3;) denotes the exterior unit normal to dD. Since / G LP{D), 
we have u G W'^’'p{D) [26| and u G by embedding [2|. Since n = 0 on dD, we have 

dnu{x + tn{x)) dt| < (e d)\\dnu\\Loo^D\D,+s-) < C{e + 6 )\\u\\w 2 ,p^d^, 
where we used |(ir(®)| < e + 5. The other assertion is shown similarly. □ 
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5. GALERKIN APPROXIMATIONS 


The error analysis of the Galerkin method defined in Q can be divided in two parts, 
i.e., similar to the previous sections, it is sufficient to show error estimates in Di and in 
D 2 separately. In the following we will again only consider the case D = Di and = 

Di \ Sf^, while the case D = D 2 can be treated similarly. As in the introduction, we let 
Vh = Pi{Th) n Hq{Q,) be the finite element space of continuous piecewise affine functions 
associated to Th- Furthermore, we denote by Ih '■ G°(n) —)• Pi{Th) H H^{^) the nodal 
interpolation operator. By construction of we see that {T G Th : T n / 0} is a shape 
regular triangulation of Df^. The corresponding finite element space is obtained by restriction, 
i.e. Wh = {vh\Di ■ Vh G Vh}. Since g G C°{Q), we can define 

V^h,g = {vh GWh:vh = hg on dDl}. (19) 

Compatible with Q , we define G ^ by 


>Di 


Vtt| • Vufe dx = / fvh dx for all Vh G qj 


'Df, 


( 20 ) 

TeTh] 


and set uf^ = IhQ on D\Df^. Furthermore, we let k = max{diam(T) : TndS^^^ 7^ 
and h = max{diam(T) : T G Th}- We then have the following statement. 

Theorem 5.1. Let G be the solution to (16) such that = g on dDf^, and let 

u\ G ^hg be the solution to (20). Then 

— tl|||i7l(_D) < C{\f5 + K3 + /i)(||/||l2(^) + \\g\\H^(D))- 

Proof. We observe that uf^ — Vh G W^q C HQ{Dff) for any Vh G Whg- Therefore, Cea’s Lemma 
m provides the following quasi-optimal approximation result 


< C inf \\u^’^ - VhWn'^iDi)- 




( 21 ) 


We will estimate the best-approximation error in the following. Let G be the 

solution to (10) with Dirichlet boundary datum g on dD^, and let G be the 

solution to (10) with replaced by and Dirichlet boundary datum g on In 

view of (Pl^we set = g on D\ , which implies G H^{D). Using Theorem |3.1 


Theorem 4.1 with f = f + Ag, we obtain 


inf 




inf 


- Vh\\H'^{Dn < 




Vh\\m{Dl) + — V^\\m{Dl) + 11'“'^ “ 


< inf - Vh\\m(Di) + CV6{\\f\\L2(^D) + \\g\\H2(D))- 


Vh&W, 




By embedding 0, we have that G C^{D). Therefore, IhU^^^ G is well-defined and 


inf 




ffi" ' “ — < ||rt^ — IhU^ \\m{Dl)- 


Next, we estimate the right-hand side of the latter estimate on each element T G Th. We 
distinguish three cases. 

(i) T n = 0. Since = g\x G H^{T), standard interpolation error estimates [T5l 
Theorem 4.4.4] yield 

- IhT^^\\H-L{T) = \\g - hglln'^iT) < C'diam(r)||tt||j^2('r). 
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(ii) T C . Then G H‘^{T), and as in (i) we have 

11< C'diam(r)||M^+'^||^2(7^). 

(hi) int{T)f^^D^~^^ / 0. By embedding [21 Theorem 5.4], we have that G 
for |?<ooifn = 2orp = 6ifn = 3. Since = g on ^ we further have that 

^£+<5 g IT^’^(T). Using Holder’s inequality and |15l Theorem 4.4.4] we thus obtain 

< \T\'^ — IhU^~^^\\y/i,p{T) < C\T\i\\u^~^^\\wi,p{T) 

^ C\T\3 {\\g\\fj2(^X\D^+/>) + \W^~^^\\H^{TnD^+S))- 

In cases (i) and (ii) we have that diam(T) < h, and in case (hi) we have that |r| < diam(T)^ < 
K^. Therefore, collecting the estimates in (i), (ii) and (ih), we obtain 




TcDf 



HMD ~b ^ ~ ^hU^^^\\\i{T) 

int{T)r\dD^+^=9 

ml(T)naD®+^^0 

<C ^ ^2(2^)-1-C 

5 i22(2’\Ds+^) “b \\H^{TnD^+^)') 

TcD^+^ , 

int(r)naD^+V0 

< C{h? -1- K3) 7i^+‘^ ^2(£)£+«) + 5 ^2p£\£)£+«)- 


Since \\u^^^\\H 2 (^r>e+s) < C{\\f\\L 2 (D) + \\g\\H2(D)) by Lemma |2^ and = I^g - g on 

D \ which can be estimated as above, the proof is complete. □ 

Remark 5 . 2 . According to [ 29 ] we have G H^/^+®°(L)|) n with sq > 0 depending 

on the shape regularity of Th- The above error estimates may alternatively be derived by 
estimating inf^gu^e^ using interpolation, cf. [T 5 |. However, the constants will 

depend on \\u'^’^\\h3 /2(^£)^'^, and as the number of re-entrant corners in Df^ might in general be 
unbounded as /i —?■ 0 the corresponding estimates might not be uniform in h anymore. 


L^-estimates are derived again via a duality argument. Note that the Dirichlet problem 
on D‘^ is not H^-regular, and thus, we cannot rely on standard estimates as in the previous 
sections. Duality arguments for the approximation of inhomogeneous Dirichlet boundary 
value problems for H^-regular problems have e.g. been investigated in [8|. 


Theorem 5.3. Let G H^(D^) be the solution to (16) 
ul G g be the solution to (20). Then 


with = g on dDf^, and let 


,£,h 


- 'u|||l2(£,) < C{S + K3 + h'^){\\f\\L'2{D) + \\9\\m(D))- 


Proof. We define — tt|, and let G Hq{D^) n be the solution to 


Then, we obtain upon integration by parts 

Jd^ JdD' 

We estimate the terms on the right-hand side separately. 


<9, 


e,h £,h 
c 


da. 


( 22 ) 
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(i) Using Galerkin orthogonality, we obtain 

[ dx + [ dx 

Jd^ Jd’^\di Jdi 

< ||V 2 ;^’^||i 2 ('£) 6 \£)£)||Ve^’^|| 2 , 2 (£)e\ 2 ?^) + ||Ve'^’^||£, 2 (£)£) inf \\Vz^'^ — Vvh\\L' 2 {Di)- 


la 


Using Theorem |A. 1 [ Lemma |A. 3 | and Lemma | 2 . 2 [ we see that 

< C'V5||e^’^||2,2(£)e). 

Dehning the set 

a = u T, 

Terh:rn95®/0 

we see that \ Df^ C Therefore, since = IhQ — 5 on \ Df^, we obtain using 
standard interpolation estimates and the definition of 5 that 

< W^Ihg — '^g\\L'^{Cl) ^ C'5||5'||j^2(£e), 


(ib) The term ||Ve^’^||^ 2 (£)£) can be estimated by Theorem 5.1 The best-approximation 
error can be estimated as in the proof of Theorem | 5 . 1 [ i.e. 


V — Vi;/i||j;^ 2(£,|) < C'-\/5||e^’^||^2(£)) + C{V6 + + h)\\e‘^'^\\]^2(^jjy 


Summarizing, for the hrst term on the right-hand side in ( 22 ) we have 

Vz"’'‘Ve"’'‘dx < C{Vd+ !<,'' 

I 

(ii) An application of the Cauchy-Schwarz inequality and LT^-regularity of z^'^ yield 


Vz^’^Ve^’^ dx < C{Vd + + ^)^(II/IIl 2 (£)) -|- \\g\\H^{D))\\^‘^’^\\L^(D)- ( 23 ) 


IdD^ 


dnZ^'^e^’^ da < ||c?„Z^’^||i 2 ( 9 £)e) ||e^’^||i 2 (g£)e) < C'||e^’^||i 2 (£)e)||e^’^||i 2 (g£)e), 


where we have also applied Lemma A .3 We now use = Ihg — g on C\ and dD'^ C to 
estimate further 


^c’'Xhsd.)= E ll''.9-9ll 


L‘2(dD'^f\T) 


TcCt 


< 


Y, |5zi^nr|||45-ff|lioc(r) 


TcCl 


< 


Y, |5zi-nr|diam(r)||5||^2(r) 


TcCt 


<Cd^\\g\\H^(ciy 

Here, we used that \dD‘^ n T| < C'diam(T) and diam(T) < 6 for all T such that T 
and Hi Corollary 4 . 4 . 7 ] in order to estimate the interpolation error. 


Summarizing, for the second term in ( 22 ) we have 


'dD^ 




( 24 ) 
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The proof is finished by observing that 


|e s ||e 


\\l'^(d\d^) ^ Ip IIl2(_d=) 


+ Ch^ 




and using (22)-(24) and Young’s inequality to estimate the first term on the right-hand side 
of the latter inequality. □ 


A combination of Theorem 1.2 and Theorems 5.1, 5.3 proves Theorem 1.3 


Remark 5.4. Instead of incorporating the condition = g on Sf^ in the approximation 
space, this condition might also be incorporated via a saddle-point formulation. For instance, 
one might consider: Find (u^,A^) G x Xf^ such that 

/ + / Xlyhdx = 

Jn Jsf^ 

/ uliih dx = 

Jsi JSI 

for all {vh, f^h) G V/i X Xf^ with Vh as above and Xf^ chosen such that an inf-sup condition 
holds, i.e. there exists c > 0 such that for all Hh G Xf^ 


fvh dx. 

(25) 

gHh dx 

(26) 


RhVhdx > c\\vh\\m{n)\\Rh\\xi for all Vh G 14 


(27) 


holds. On the continuous level such a condition may be verified for functions /r G Y being 
defined as the closure of with respect to the norm 


x= sup / g.vdx/\\v\\Hi{Q)- 


The verification of (27) is completed by constructing a Fortin projector P : —)■ Vh 

such that P is bounded and 

/ g,h{Pv — u) dx = 0 for all v G H^{Sh) and ^h £ ^h, 


see m Proposition 2.8]. If Xh = • ^h G 14}, the construction of a Fortin projector 

amounts to 77^-stability of the L^-projection, which has recently been shown in [6] for a large 
class of finite element approximation spaces. We remark that, if g is replaced by Ihg in (26), 


then the choice /r/i = — Ihg yields n| = Ihg- Choosing u/j G Vj^Q then shows that n| is a 

solution to (j^. Penalized saddle-point problems might also be considered, cf. |16[ II.4], and, 
for instance, the method in [l3] might be interpreted as such. Moreover, surface PDFs, when 
appropriately extended to S^, might be incorporated as an additional constraint. 

6. Numerical Examples 

Demonstrating the validity of the above derived results, we set D = (—2, 2) x (—2, 2) C 
and let P = 5i3i(0). Furthermore, we let x = (xi, X2) G D and 

u{x) = {A - xl){A - xI)xd 2 {x) + {^ - xl){A - xl) - |xp)x;^(x), 

g{x) = (4 - Xi)(4 - xl) cos(l - jxp), 

/(x) = — Au(x), X G D \ F, /(x) = 0, x G F. 
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Figure 2. Sketch of computational domain and corresponding triangulation. 
The solid black line represents T. The dotted black lines correspond to dS^. 
The gray area corresponds to Sf^ for e = 1/8. Left: h = 5 and 361 vertices in 
total. Right: 5 = h? and 593 vertices in total. 


Here, Xd^ arid xd 2 denote the characteristic functions of Di and D 2 , respectively. One 
verifies that this choice of functions yields a solution to 0 . For the sake of simplicity, 
we have chosen an arbitrary globally defined and sufficiently smooth function g such that 
g{x) = u{x) for x G F. As noticed in the introduction, if g is defined on F only, we have 
to construct a suitable extension of 5 ^ to a neighborhood of F. For our numerical tests, we 
employed the linear nodal interpolant Itdr of (ir(x) = |x| — 1 in order to define 5^, namely, 
setting w|(x) = ^(1 + 5'(—we note that supp|Va;|| = Sf^. In all our tests below, we 
do not employ aligned meshes, i.e. the interface approximations dSf^ are not aligned with F. 
In Figure we have depicted the geometric setup. The Galerkin approximation uf^ is then 
computed via Q with / replaced by Ihf, which is a sufficient approximation for / on II \ 5^ 
for the chosen /. 

In our first experiment we choose a uniform triangulation, i.e. 5 = h = k, with 332 929 
vertices; cf. Figure]^ left. The convergence results for different values of e = 1/2*, i = 1,..., 20 
are depicted in Figure We observe the predicted convergence rates 0{e) for the L^(n)- 
norm, p G {2, 00 }, and O(x^) for LI^(n)-norm. In particular, the error behaves monotonically 
and saturates at a certain level, which is due to the chosen mesh. 

In a second experiment we also chose uniform triangulations, i.e. 6 = h = k, and e = 1/2^^ 
fixed. We used different mesh sizes with number of vertices in {5 329, 21025, 83 521,332 929}. 
Notice that in this example h > 10“^ and e ~ 10“®, i.e. s ^ h. The convergence results 
for different mesh sizes are depicted in Figure [^We observe the predicted convergence rates 
0{h) for the LP(n)-norm, p G { 2 , 00 }, and 0{-vh) for the LI^(II)-norm. 

In a third experiment we chose a locally refined mesh such that 5 = h?‘ and k < 46, 
which has been obtained from the meshes in the second experiment by repeated refinement 
of those triangles having nonempty intersection with dS^-, cf. Figure right. The resulting 
meshes have number of vertices in {10 053, 41445,168 717,681957}. The diffuse interface 
width e = 1/2^^ is fixed. The convergence results for different mesh sizes are depicted in 
Figure]^ We observe the predicted convergence rates 0{h‘^) for the LP(Il)-norm, p G (2, 00 }, 
and 0{h) for the 77^(II)-norm. 
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Figure 3. Convergence rates in s: Uniform mesh with 332 929 vertices. The 
solid lines correspond to the predicted rates (e, ^/£, s from left to right). The 
actual errors ||m - \\u - ul\\, Ik - (from left to right) 

are marked by crosses. 





Figure 4. Convergence rates in h: Uniform mesh with number of vertices in 
{5 329, 21 025, 83 521, 332 929} and fixed s = 1/2^*^. The solid lines correspond 
to the predicted rates {h, y/h, h from left to right). The actual errors |k — 
'«|IIl2(u)> \\u - ulW hi( n), Ik - <||l“(u) (from left to right) are marked by 
crosses. 


All these convergence rates are predicted by Theorem |1.3 
estimates on the continuous level. 


cf. Theorem 4.3 for the L°°- 


For locally refined meshes around dS^ such that <5 = 
we thus recover the convergence rates of the usual finite element method when used in 
combination with aligned grids. The example presented here also shows that these convergence 
rates are sharp in general. The number of degrees of freedom roughly doubles in two spatial 
dimensions compared to the corresponding uniform meshes, i.e. the computational complexity 
increases only by a multiplicative factor independent of the degrees of freedom. For n = 3, 
the situation is worse regarding the number of degrees of freedom and the results presented 
here might be extended using anisotropic finite elements, see e.g. [3]. 


7. Conclusions 

In this paper a diffuse interface method for solving Poisson’s equation with embedded in¬ 
terface conditions has been investigated. The diffuse interface method could be interpreted 
as a standard Dirichlet problem on approximating domains. Thus, this method can also be 
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Figure 5. Convergence rates in h: Locally refined mesh with number of 
vertices in {10 053,41445,168 717,681957} and fixed e = 1/2^*^. The solid 
lines correspond to the predicted rates (h^, h, h? from left to right). The 
actual errors ||u - \\u - uI\\h^q), \\u - u|||x,oo(j^) (from left to right) 

are marked by crosses. 


used to numerically solve Poisson’s equation with Dirichlet boundary conditions on compli¬ 
cated domains. Error estimates in L?- and L°°-norms have been proven and verihed 
numerically. The use of uniform meshes gave suboptimal convergence rates in terms of the 
mesh size. This is explained by the fact that the approximating domains are polygonal and 
non-convex in general, and therefore the corresponding Dirichlet problems do not allow for 
Lf^-regular solutions in general. Using locally refined meshes we could recover order-optimal 
convergence rates in terms of the mesh size. In two spatial dimensions, the computational 
complexity is only increased by a constant factor when using locally refined meshes. For three 
dimensional problems, the analysis might be extended to cover anisotropic finite elements, 
which leads to a computationally efficient method. The presented method might furthermore 
be applied to more general second order elliptic, and the generalization to parabolic equations 
with r = r(t) seems possible. Moreover, the method might be extended to interface problems 
where g = g{u) couples the two subproblems investigated here. It is open, whether one can 
improve the method, for instance by a post-processing step, in order to retrieve order-optimal 
convergence with respect to the mesh size (at least away from the interface) in the case when 
quasi-uniform grids are used. 
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Appendix A. Basic Properties of Diffuse Approximations 


We let D = Di and P = dD in the following. The corresponding estimates for D 2 are 
derived in a similar fashion. Let Sq > 0, 0 < e < eo and 


= {x C D : dist(x, P) > e}, Pg = D \ . 

For £0 sufficiently small, the projection pp : Fgp —)■ P is well-defined, and given by [22l Chapter 
7, Theorem 3.1] 

pr{x) = X — dT{x)Vdr{x). (28) 


Notice that Vdr(T) = n{pQp){x)) for all x G Tg^ [221 Chapter 7, Theorem 8.5], and we will 
therefore write n{x) = Vdr{x). Here, n denotes the exterior unit normal field to dD. For 
t G [0,eo], we define the mapping d'i(x) = x — tn{x), x G F, and note that ‘^^(F) = dD^. 
Moreover, D<ht(x) = I — tD'^drix), and, cf. [I?! Eq. (9)], 

sup j det D^t{x) — 1 + tA(ip)(x)j < Ct^. (29) 

xSF 

Hence, after decreasing ep if necessary, we may assume that 

- < det <2, 0 < t < Eq. (30) 


For any integrable v the transformation formula implies 



v{x) dx 


v{x 


tn{x))\ det D^t{x)\ dt dcr(x). 


r Jo 


(31) 
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For the right-hand side we will employ the fundamental theorem of calculus 
v{x — tn{x)) = v{x) — / dnv{x — sn{x)) ds, 0 < t < eo- 

Jo 

The following is a central estimate; cf. m Theorem A.2], 

Theorem A.l. There exists a constant C > 0 not depending on e such that 
||^i||L2(r,) < C'(\/e|kllL2(r) + e||'9nf ||L2(r,)) for v G H\Te). 
Proof. Let 0 < t < e. Using (|32|) and Holder’s inequality we obtain 


(32) 


\v{x — tn{x))\‘^ < 2{\v{x)\‘^ + e / — sn(x))p ds). 

Jo 


Using the latter in (31) and using (30), we obtain 


l^llL2(r^) < C'(£||^llL2(r) + e \dnv{x — sn(x))p dsdfj). 


Using (31) and (30) we further may write 

[ [ \dnv{x — sn{x))\^ ds da < C f \dnv{x)\‘^ dx, 

Jv Jo Jv^ 

which completes the proof. □ 

Lemma A. 2. There exists a constant C > 0 independent of e such that for v G 

\\'o\\L^(dD^) < C'(ll'f^llL2(r) + Vf\\dnV\\L-2(T^))- 
Proof. The transformation formula yields 

/ |upd(T= / |u(x — en(x))|^| det L><h£(x)| d(T(x). 

JdD^ Jr 

Using (32), the assertion follows with similar arguments as in the proof of Theorem A.l □ 
Lemma A. 3. There is a constant C > 0 independent of e such that for v G 

\\'0\\L^{dD‘^) < C'||U||^1(£)£:). 

Proof. We extend v G H^{D^) to a function v G by reflection [H Theorem 4.26]. We 

have ||'c||j 71 (r™) < C\\v\\Hi(^Dsy Continuity of the embedding H^{D) ^ Lp‘{dD) [2] implies 
ll'^llL 2 (ai 5 ) < C||u||j:^i(£)). Lemma A.2 and Ur'IlHqrg) < complete the proof. □ 

The preceding lemma with a different proof may also be found in [U Lemma 2.1]. 









